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ABSTRACT 

Context. Most semi-empirical static ID models of the solar atmosphere in the magnetically quiet Sun (QS) predict a rise of the tem- 
perature at chromospheric layers. Numerical simulations of the solar chromosphere with a variable degree of sophistication, i.e., from 
one-dimensional (ID) to 3D simulations, and assuming local thermal equilibrium (LTE) or non-LTE (NLTE), on the other hand only 
yielded an increase of the brightness temperature without any stationary increase of the gas temperature. 

Aims. We investigate the thermal structure in the solar chromosphere as derived from an LTE inversion of observed Ca n H spectra in 
QS and active regions (ARs). 

Methods. We applied an inversion strategy based on the SIR code to Can H spectra to obtain ID temperature stratifications. We 
investigate the obtained temperature stratifications on differences between magnetic and field-free regions in the QS, and on differ- 
ences between QS and ARs. We determine the energy content of individual calcium bright grains (BGs) as specific candidates of 
chromospheric heating events. We compare observed with synthetic NLTE spectra to estimate the significance of the LTE inversion 
results. 

Results. The fluctuations of observed intensities yield a variable temperature structure with spatio-temporal rms fluctuations below 
100 K in the photosphere and 200-300 K in the QS chromosphere. The average temperature stratification in the QS does not exhibit 
a clear chromospheric temperature rise, opposite to the AR case. We find a characteristic energy content of about 7 x 10 18 J for BGs 
that repeat with a cadence of about 160 sees. The precursors of BGs have a vertical extent of about 200 km and a horizontal extent 
of about 1 Mm. The direct comparison of observed with synthetic NLTE profiles partly confirms the results of the LTE inversion that 
the solar chromosphere in the QS oscillates between an atmosphere in radiative equilibrium and one with a moderate chromospheric 
temperature rise. Two-dimensional x - z temperature maps exhibit nearly horizontal canopy-like structures with a few Mm extent 
around photospheric magnetic field concentrations at a height of about 600 km. 

Conclusions. The large difference between QS regions and ARs, and the better match of AR and non-LTE reference spectra suggest 
that magnetic heating processes are more important than commonly assumed. The temperature fluctuations in QS derived by the LTE 
inversion do not suffice on average to maintain a stationary chromospheric temperature rise. The spatially and vertically resolved 
information on the temperature structure allows one to investigate in detail the topology and evolution of the thermal structure in the 
lower solar atmosphere. 
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1 . Introduction [ Vernazza etaDll976l: lUlmschneider et al][l978l iLeenaarts etail 

201 lh . Such ID models with a temperature reversal faithfully re- 

The thermal structure of the solar chromosphere remains enig- produce observed (average) spectra to the most minute level, but 

matic and to the greatest extent unexplained more than a century a ll lack something: they provide a temperature stratification with 

jJ ■ after t he first descr i ption of the ex istence of the chromosphere height in the solar atmosphere, but do not explain the reason for 

itself (ISecchill 18601: ILockyeill 18681) . Nowadays, two competing its shape. 

views on the thermal structure exist, i.e., a static/stationary and Trie dynamical view of the solar chromosphere is com- 
a dynamical view. parably newer than the static ID models and resulted thanks 
The static view is expressed by a long series of one- to the increase in the available computing power during the 
dimensional (ID) t hermal stratifications that reproduce ob- last two decades. Opposite to the static models that were de- 
serve d spectra (e.g.. iGingerich et al.lll97lt iHolweger & Mulled rived from observations, the dynamical approach usually goes 
19741: IVernazza et al.l 119811: iFontenla et al.l Il993l 120061 see the opposite way and tries to model the solar chromosphere 



Rutten & UitenbroekT (1201 2l) and iRuttenl (12012 ) for additional from basic principles, i.e., numerical simulations that encom- 



references). To produce emission in the line cores of chromo- pass the differential equations that gove rn the temporal e volu- 

spheric spectral lines, all these ID models require a reversal of tion of material in the solar atmosphere (iNordlundll 19821) . One 

the temperature gradient with height at some point. Their com- of the advantages of such a forward modeling is that one only 

mon feature therefore is the existence of a temperature mini- needs to decide on what are the relevant physical processes that 

mum at a height of about 500 km above the photosphere, where must be included in the set of differential equations. Whereas 

the gas temperature starts to rise towards higher layers (e.g., numerical simulations of the solar photosphere abound (e.g., 
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Nordlund & Steinl 1 19901: [Sterner et al.l [19981 [Stein & N ordlun 



1998[ iSteiner et al.l 120081: iFabbian et al.l 120101: iBeeck et al 
2012b. there are relatively f ew ch r omospheric simu l ations 
(iRammacher & Ulmschneiderl 119921: ICarlsson & Steinl Il997b 



Wed emever et al.l 120041: ISchaffenberger et al.l 120061: ICuntz et al I 
2007t iLeenaarts et al. 200|). 



One reason is the spatial, and especially the vertical extent 
of the simulation box that has to reach up to a height of about 
2 Mm with a sufficiently fine sampling on the order of a few 
tens of km to encompass the solar chromosphere. A second 
limitation is caused by the physical conditions in the upper solar 
atmosphere. The gas density in the solar photosphere is high 
enough that the assumption of local thermal equilibrium (LTE) 
is justified. In LTE, the available total energy is distributed 
over all available degrees of freedom by frequent collisions 
between particles such that each degree of freedom obtains 
a fully determined fraction of the total energy. The exact 
distribution solely depends on the physical laws holding in 
the material of interest. Assuming LTE, the values of several 
physical quantities such as population states, thermal motions, 
etc. are determined instantly as soon as a temperature value is 
given. Because of the decrease of the gas density with height 
in the gravitationally stratified solar atmosphere, the LTE 
condition breaks down when a sufficient number of collisions 
takes longer than the dynamical evolution of the physical 
quantities that can change by, e.g., radiative proce sses (e.g., 
ICarlsson! 120071: IWedemever-Bohm & Carlssonl 1201 ll) . A com- 
plete treatment of this regime in a non-LTE (NLTE) formalism 
still exceeds the computing power of even the most recent 



computer clusters ( Auer & Mihalas 


1969; Avres & Wiedemann 


1989 


: IRammacher & Cuntz 


1991 


; Carlsson & Stein 


1992, 


1997 


; [Uitenbroek 2000 


: Asensio Ramos et al.l 


2003; 


Leenaarts et al.l 120071: IWedemever-Bohm & Carlssonl 
Carlsson & Leenaartsll2012HMartmez-Svkora et al.ll2012l). 


2011; 





The dynamical numerical simulations of the solar chromo- 
sphere revealed a quite distinct thermal structure than the static 
ID models. The brightness temperature at chromospheric lay- 
ers tended to fluctuate strongly, leading to intermittent emis- 
sion in the synthesized chromospheric spectra. The gas tem- 
perature, however, did not show any temperature reversal on a 
temporal average, contrary to most of the stat i c ID models (cf. 
ICarlsson & Stei n 1995, 1997; Kalkofen 200ll l2012l) . Thus, the 



average brightness and gas temperatures in such simulations dif- 
fered significantly, while the average synthesized chromospheric 
line spectra corresponded to neither the average gas nor the aver- 
age brightness temperature. This puzzling behaviour even raised 
the questio n how to define a "mean" temperature at all under this 
condi tions (IRammacher & C untz 20 051: lUitenbroek & Criscuolil 
l201ll) . This point is actually more important than it seems at a 
first glance because the static ID models are based on the re- 
lation between brightness temperatures (i.e., observed spectra) 
and gas temperature (i.e., model output). Questioning this rela- 
tion between gas and brightness temperature, or the definition 
of temperature therefore converts directly into questioning the 
results of the static ID models that prov ide an excellent fit to 
observed average spectra (lKalkofenll2012l) . 

There are some common problems for both the static and the 
dynamical view of the solar chromosphere. One is the width of 
spectral lines synthesized from the simulations or models, re- 
spectively. Observed chromospheric spectral lines are usually 
found to be significantly broader than the synthetic lines, (far) 
in excess of the thermal line width for any reasonable tempera- 
ture range. This behavio ur shows up most prominently in off - 
limb observations (e.g.. iBeckersI 119681: iBeck & Rezaeill201ll) . 



Part of this excess in line width can be explained by the spa- 
tial resolution of the observations that represent a spatial aver- 
age over solar structures with different physical properties, and 
especially different macroscopic mass flows with their related 
Doppler shifts. But also the average of synthetic spectra from 
simulations with spatial grid sizes down to a few ten kms lacks 
the observed line width, which indicates that some physical in- 
gredient is missing in the simulations, or that their spatial sam- 
pling - despite being below the photon mean free-path length - 
still is insufficient to capture all of the dynamical evolution of, 
e.g., small-scale mass flows. The line width in both static and 
dynamical models is therefore often adjusted ad-hoc by some 
additional broadenin g parameters (cf . lAUende Prieto et al1l2 004: 
lAsplund et al .1120041 and references therein), usually labeled as 
macroturbulent and microturbulent velocities that represent spa- 
tially un resolved line-broadening m echanisms of unknown ori- 
gin fe.g. JCrad[T^lRuttenll200l . 

A second problem is how to reconcile the temperature strat- 
ifications with t he observations of molecu lar lines at chromo- 
spheric heights (lAyres & Testermanl Il98ll) . This problem af- 
fects more the static than the dynamical models, because in 
the former the location of the temperature reversal at low 
heights prev ents any significant formation of molecules at l ay- 
ers above it (lAvresiri98lb iKalkofen et aDll999k lAvresll2002h . A 
"cool" chromospheric part could partly be included also in static 
models by "post poning" the chromospheric temperature rise 
to higher layers (IAvresll2002l lAvres etal.ll2006l iFontenla etaD 
120071: iRezaei et al . 2008), but it was imposed in some sense as 
a post-facto crutch. In most numerical simulations, phases with 
a temperature that is sufficiently low to allo w the formation of 
molecules were usually auto matically present (I Wedemever et al.l 
120041: ILeenaarts et al.ll20TTb . 

A third problem is a discrepancy between observations in 
different spectral lines that yield seemingly contradictory re- 
sults. Observations with the Solar Ultraviolet Measurement of 
Emitted Radiation (SUMER) instrument on-board the SOlar and 
Heliospheric Observatory (SOHO) satellite in ultraviolet (UV) 
lines revealed that all these lines are in permanent and ubiquitous 
emission wit h little intrinsic variati on in quiet Sun regions of the 
solar surface (ICarlsson et al.ll 19971) . All these lines are classified 
as "chromospheric". Contrary to this, chromospheric Ca II lines, 
especially H and K in the near-UV, show large spatio-temporal 
fluctuations from reversal-f ree absorption profi l es to moderate 
emission in the line cores (|Liu & Smith] Il972t ICram & Darnel 
119831: iRezaei et al.ll2008l: IBeck e t al. 200j|. While the temporal 
evolution is naturally explained in dynamical simulations, per- 
sistent emission is only reproduced in static ID models with a 
permanent chromospheric temperature rise. 

The way to resolve this contradictory behaviour of observed 
spectral lines and to reconcile the conflicting views of the 
static and dynamical solar chromospheric models is not clear 
at present. There is one additional information to be considered: 
the temperature in the solar photosphere is about 6000 K, while 
the ion and electron temperature in the solar corona exceeds one 
million K, where both results are undoubted. Thus, there has to 
exist a reversal of the temperature gradient at some height, but 
where exactly is unclear. 

Magnetic fields could be a possible mediator between 
the two conflicting scenarios of the solar chromosphere. 
Firstly, magnetic fields have not been considered neither in 
all of the static nor many of the dynamical models (but se e 
also iMartmez-Svkora et al.ll2012b iKhomenko & Colladosll2012l) . 
Secondly, in quiet Sun (QS) conditions, concentrations of mag- 
netic flux are theoretically expected to spread out from their con- 
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fined diameter in the low photosphere to fill the full volume in 
the chromosphere. This structure of several expanding flux tubes 
creates automatically a dividing layer at some chromospheric 
height by the so-called ma gnetic canopy tha t separates field-free 
from magnetized plasma (iGiova nelli 1980; Stein er et"aL| [l986: 
Solanki et al] 1 199 It IStepan & Truiillo B uenol 120101: iNutto et all 



201J Holzreuter & Solanki 2012b . i.e., that separates a dynam 
ical chromospheric regime below from a stationary one above 
the canopy (as sketched for instance in iBeck et aD 120081 their 
Fig. 22). 

In this contribution, we investigate the spatio-temporal be- 
haviour of the thermal structure in the solar chromosphere as 
derived from an inversion of Can H spectra assuming LTE 
(IBeck et al.ll2013l BE13). This study follows the approach of the 
static ID models by deriving a gas temperature from an observed 
brightness temperature, but uses spatially resolved spectra that 
trace the dynamical properties of the solar chromos phere - as 
in the dynamical view of the chromosphere (see also iHenriquesI 
l2012h . We use observations of both QS and active regions (ARs) 
to investigate differences between these types of solar structures. 
The observations and the inversion scheme are described in de- 
tail in BE13, we thus refer the reader to the latter publication. 
We will only highlight some observational limitations that are 
relevant for the current contribution in Sect. [2] We then concen- 
trate on the thermal structure of the chromosphere in QS and 
ARs in Sect. [3 both on average and spatially resolved. The find- 
ings are discussed in Sect.|4] Section[5]provides our conclusions. 
Appendix lAl discusses the possible amount of spectral broaden- 
ing by different mechanisms acting on the observed spectra. 

2. Relevant observational details 

We use in total five different observation s taken with th e 
POlarimetric LIttrow Spectrograph (POLIS, IBeck et al.ll2005h . 
POLIS deliverecO intensity spectra in Can H at 396.85 nm to- 
gether with spectropolarimetric data of the two Fei lines at 
630 nm. Two data sets were taken in QS regions, i.e., a time 
series of about one-hou r duration and a large-area map (cf. 
IBeck et al.1 12009L I2012L BE09 and BE12 in the following). 
The other three data s ets are large-area maps of the active re- 
gion NOAA 10978 fcf. lBeck & Rammache r feoiOtlBethge et al.1 
120121) . The parameters most relevant in the present context are 
the slit width of 075, the step width in the case of spatial scan- 
ning of also 075, the temporal cadence of the time series of 
21 sees, and the spatial sampling along the slit of 07292 in the 
Can H spectra. These values provide a pixel scale of 075 x 
07292 in the Ca spectra. The data were taken with real-time see- 
ing c orrection by the Kiepenh euer-Institut adaptive optics sys- 
tem (Ivon der Luhe et al.1 120031) and have a seeing-limited spa- 
tial resolution of about 1", which roughly coincides with the 
sampling-limited resolution in the scanning direction. 

Besides from the study of individual events in the time series 
(Sect. 13.5b , always the average of the statistics in all individual 
maps is used to expand the statistical base. 



3. Results of LTE inversion 

The inversion setup for Ca 11 H described in BE 13 yields temper- 
ature stratifications for spatially resolved spectra. With the inver- 
sion results, we now can investigate the statistics of the temper- 
ature distributions at different optical depth levels log r, analo- 
gously to BE 12, where the intensity distributions at different X 




POLIS was officially de-commissioned end of 2010. 
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Fig. 1. Relative occurrence of temperatures in the QS. Left top: 
temperature versus (vs.) optical depth for the full FOV. The rela- 
tive occurrence of a given temperature is indicated by the colour 
with the corresponding colour bar at top. Thick black line: orig- 
inal HSRA model. Dashed', average temperature. The vertical 
dashed lines denote the layers of log r for which the temperature 
histograms are shown in the rest of the panels. In the histograms, 
black/red/blue lines denote the full FOV/magnetic/field-free lo- 
cations. 



and the characteristic properties of different spatial samples in 
the observed field of view (FOV), i.e., field-free or magnetic re- 
gions, and the full FOV, were analyzed. The distinction between 
field-free and magnetic sample is made using the simultaneous 
polarimetric observations in the 630 nm channel of POLIS by 
setting a threshold in the observed polarization signal (cf. BE12). 

3.1. Average temperature and temperature statistics in QS 

Following the approach of BE 12, we calculated the characteris- 
tics of the temperature distributions up to second order (mean, 
root-mean- squared (rms) fluctuation, skewness). Figure Q] shows 
some of the statistics of temperature versus (vs.) optical depth 
for the maps taken in the QS. The upper left panel shows the 
relative occurrence (denoted by the colour bar) of temperatures 
in the full FOV for each optical depth level. The three sam- 
ples in the FOV (field-free, magnetic, full FOV) differed only 
slightly in this way of presentation, therefore only the latter 
is shown. The remaining panels show the histograms of tem- 
perature for all spatial samples at three optical depth levels of 
logr = -0.6, -2.6, and -4.6, respectively. The magnetic loca- 
tions (blue lines) show a weak preference for slightly higher 
temperatures compared to the full FOV or the field-free loca- 
tions. The average temperatures of all three samples do not show 
an indication of a chromospheric temperature rise as predicted, 
e.g., by the orig inal Harvard Smithsoni an Reference Atmosphere 
model (HSRA. iGingerich et al]|l97lh . e.g., compare the black 
solid line (HSRA) and black dotted line (< T > in the full FOV) 
in the upper left panel of Fig. [T] 

The lower panel of Fig. [2] shows the variation of the rms 
fluctuations of temperature with optical depth and geometrical 
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Fig. 2. Bottom panel: rms fluctuations of temperature for the full 
FOV (black), field-free locations (red), and magnetic locations 
(blue). The three dotted lines indicate the maximum deviation 
from the average temperature. The dashed line denotes the rms 
fluctuation in an AR map. Top panel: skewness of the tempera- 
ture distributions. The vertical dotted line indicates z = 150 km. 



height. For the conversion from one to the other we used the 
tabulated relation given in the HSRA atmosphere model, oppo- 
site to a derivatio n of the height scale from the in version results 
themselves as in PuschmannetaD (2005, 20lo|). We note that 
in the derivation of the temperature statistics always the indi- 
vidual temperature stratifications vs. optical depth were used, 
and that the rms value will (slightly) underestimate the char- 
acteristic range of fluctuations in case of asymmetric distribu- 
tions (cf. Fig. []]). The geometrical heights in Fig. [2] only in- 
dicate a probable value for the average location of the forma- 
tion heights. As in the analysis of the intensity statistics of Can 
H and Can IR spectra in BE 12, we find that the rms fluctua- 
tions increase nearly monotonically in geometrical height, with 
a minimum of fluctuations near z ~ 100 km (corresponding to 

log r 1) caused by t he temperature inversion of gr anules and 

inter-granular lanes (cf . iPuschmann et al1l2003ll2005l) . The fluc- 
tuations stay below 100 K up to 500 km height, and reach only 
200- 300 K at 1000 km. The magnetic regions show slightly 
higher rms fluctuations at all heights and pass through the mini- 
mum of rms fluctuations lower in the atmosphere than the other 
two samples. The latter behaviour is presumably caused by the 
shift of the optical depth scale on magnetic locations. The max- 
imal temperature variations (dotted lines in the lower panel of 
Fig. [2]) reach about 1000 K at logr ~ -4 (z ~ 700 km). 

BE 12 found a pronounced difference in the skewness of in- 
tensity distributions between field-free and magnetic locations, 
i.e., an increased skewness on magnetic locations at intermediate 
wavelengths between the line wing and the core, corresponding 
to formation heights of about z =135-150 km. The upper panel 
of Fig. [2] shows that this also happens at about the same height in 
the case of the temperature retrieved by the inversion. A positive 
skewness indicates that temperature excursions towards values 
above the average are more frequent than those towards values 
below it. The skewness on magnetic locations increases strongly 
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Fig. 3. Relative occurrence of temperatures in an AR. Dash- 
dotted: average temperature in the AR. Dotted: average temper- 
ature in QS. Thick black line: original HSRA model. 



already at lower optical depths of about log r 4 when com- 
pared with the other two samples. This could again be due to a 
shift of the optical depth scale on magnetic locations, but could 
also indicate that propagating waves turn into shock fronts at 
lower heights in the rarefied medium of magnetic flux concen- 
trations than in the field-free surroundings. 

3.2. Average temperature and temperature statistics in AR 

One could argue that with the list of assumptions used in the 
inversion (LTE, unmodified HSRA density, lack of hydrostatic 
equilibrium, complete frequency redistribution), it would be im- 
possible to obtain a chromospheric temperature rise in the inver- 
sion of spectra because the (partial) decoupling between temper- 
ature and emitted intensity that exists in NLTE would be needed. 
Figure [3] shows that this is not the case, but that a chromospheric 
temperature rise results also in LTE if input spectra different to 
those in QS are analyzed. Both the average temperature and the 
distribution of temperatures in the AR data cover the temper- 
ature rise in the original HSRA model (Fig. [3j. The offset of 
the temperature stratifications in the optical depth scale - the 
HSRA temperature starts to rise at a larger height (log r = -4) 
than the average AR temperature retrieved from the inversion 

(logr 3) - is presumably caused by the off-centre location 

of the AR (6 = 50°) and the related shift of the optical depth 
scale because of the inclined line of sight (LOS). We remark 
that the drop of temperature at the uppermost optical depth lev- 
els (logr < -5.5) is purely artificial. It results from the need of 
having a reduction of temperature at the uppermost atmosphere 
layers for generating double reversals in the Ca line core when 
using the LTE approximation. The rms fluctuation of the tem- 
perature in one of the AR maps is over-plotted in the bottom 
panel of Fig. [2] It reflects the intrinsically larger spatial variation 
of temperature in an AR compared to that in QS, but we note 
that the range goes up to beyond 1000 K rms, implying even 
larger maximal variations. Temperature fluctuations of this size 
are large enough to rise the temperature from an atmosphere in 
radiative equilibrium with about 4000 K in the chromosphere to 
the temperature values of commonly used atmospheric models 
such as the HSRA. 

3.3. Correlation between intensity and temperature 

In our inversion approach, each wavelength A is attributed a cor- 
responding value in r, i.e., to adjust the intensity 1(A), the cor- 
responding temperature T(r) should be varied. The iterative in- 
version process, however, does not only modify the temperature 
at the given value of r, but uses a fourth order polynomial for 
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Fig. 4. Correlation coefficient between intensity at a given wave- 
length 1(A) and temperature r(logr). 



modifying the full temperature stratification. In addition, when 
the radiative transfer equation through the model atmosphere 
is solved in the spectral synthesis, the temperature values of 
some extended optical depth range contribute to the intensity 
at one given wavelength. Therefore, the one-to-one correlation 
between a given wavelength and its attributed optical depth is 
partially broken up. Additionally, some solar surface structures 
such as (vertical) magnetic fields leave a signature in intensities, 
and hence temperature over some extended height/optical depth 
range. 

To determine which layer in log r corresponds to the inten- 
sity at a given w avelength, we calculated th e linear correlation 
coefficient C (cf. iBeck & Rammacher1l2010l) between intensity 
vs. wavelength and temperature vs. log r. Figure |4] shows the re- 
sult for the large-area scan on disc centre. The correlation coeffi- 
cient between intensity and temperature is similar to the intensity 
contribution function (e.g., Fig. 1 of BE09) in the line wing up to 
about 396.65 nm, with a single maximum of correlation between 
log r « -0.5 to - 1 .5 that slowly moves to upper layers for wave- 
lengths closer to the line core. Between 396.65 nm and 396.8 nm, 
the correlation coefficient shows a second, slightly weaker local 
maximum of correlation at about logr « -3.5. In the line-core 
region from 396.81 nm to 396.89 nm, two clear local maxima of 
correlation exist, where the upper one above log r « -4 corre- 
sponds to the intensity contribution function, whereas the origin 
of the lower one at log r « -1 is not clear at once. 

Figure [5] reveals that the maximum of correlation between 
line-core intensity and temperature at low atmosphere layers is 
most likely caused by locations with large photospheric mag- 
netic flux. A comparison between the line-core intensity / core 
(lower left panel in Fig. 0) as the most prominent tracer of pho- 
tospheric magnetic flux and the temperature at optical depths of 
logr = -0.9, -2.5, and -5.4 shows that the photospheric mag- 
netic network can be identified at both logT = -0.9 and -5.4 
(correlation coefficient with / core of 0.50 and 0.85, respectively), 
but is completely inconspicuous at logr = -2.5 (C = 0.08). 
This implies that vertical or only weakly inclined magnetic fields 
- characteristic of the magnetic network and isolated magnetic 
flux concentrations because of the buoyancy force - leave promi- 
nent signatures in temperature in at least two height ranges, i.e., 
at upper atmospheric layers above log r = -4, and at lower at- 
mospheric layers at about log r « - 1 . This is thus most likely 
the reason for the increased skewness in the magnetic sample 
at logr « -1.2 (z ~ 150km) found in the previous section 
and in BE 12: magnetic flux concentrations and their immedi- 
ate surroundings show up brighter than the quiet inter-network 
(IN) at intermediate wavelengths, and hence have a slightly en- 
hanced temperature at the corresponding a tmospheric layers (cf. 
ISheminova et al.ll2Q05l : ISheminovall2012l) . This, however, does 
not explain why the intensity (temperature) is enhanced only at 
these wavelengths (heights), but not in between. 
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Fig. 5. Comparison of line-core intensity (lower left) and temper- 
ature at logr = -0.9 (lower right), -2.5 (upper left), and -5.4 
(upper right). The correlation coefficient between 7 core and T(r) 
is denoted in the bottom middle in the temperature panels. 

3.4. Energy flux 

3.4.1 . Temperature rms fluctuations with height 

To estimate the vertical energy flux, we determined the evolu- 
tion of the energy content with height assuming that the energy 
balance is effected in the vertical direction in a gravitationally 
stratified atmosphere with radial energy transport at least up to 
the photosphere with its convective granulation. To separate the 
dynamical part of the energy balance from any possible static 
background caused by radiative equilibrium with the incoming 
photospheric continuum radiation, we converted the rms fluctua- 
tions of the temperature Ar rms (r) - instead of using the modulus 
of the temperature T(r) - to the equivalent increase of internal 
energy using 



RAT 

Mr-D 



(i) 



where/? = 8.31 Jmol" K -1 , ji = 1.3gmol~ , and y = 5/3 are 
the gas constant, the specific mass, and the adiabatic coefficient, 
respectively (cf. Eq. (4) of BE09). The tabulated values of the 
gas density p gas in the HSRA atmosphere model were used. With 
the relation between r and geometrical height z tabulated in the 
HSRA model, this yields the (stationary) average variation of 
energy with height (top panel of Fig. [6]) . 

The resulting volume energy density E/V can be compared 
with the result of the first two papers of this series. The energy 
density derived from the direct conversion of intensity fluctua- 
tions to temperature (BE 12) lies slightly above the one derived 
from the temperature fluctuations in the inversion. The reason 
is that in the inversion the variation of intensity is implicitly 
converted to a temperature fluctuation over some optical depth 
range, whereas in the direct conversion from intensity to tem- 
perature all variation is attributed to a single optical depth point. 
The latter approach therefore leads to larger rms fluctuations 
than in the case of considering the radiative transport through 
the atmosphere. The energy density corresponding to the veloc- 
ity rms (BE09) lies below that derived from intensity fluctua- 
tions by direct conversion or by inversion up to z ~ 1 Mm, but 
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Fig. 6. Variation of energy density with height. Top: energy den- 
sity equivalents corresponding to velocity (red line with pluses), 
intensity (purple dash-dotted), and temperature rms fluctuations 
(black solid). The black dashed line results for the tempera- 
ture rms fluctuations when using the gas density in the original 
HSRA model. The blue line indicates the gas density in the mod- 
ified HSRA model. Middle: energy gradient per km in height as- 
suming a propagation with v ac . The red line results from using 
a net cooling rate of 100 W g" 1 (CL12). Bottom: radiative losses 
in the units used in CL12. 

exceeds the others above that height. We note that all temper- 
ature fluctuations - when using AT without applying Eq. (Q} - 
from either velocity, intensity, or the LTE inversion fall short of 
the energy required to lift an atmosphere in radiative equilib- 
rium to the HSRA or any similar atmospheric model (~ 1000 K 
needed at z ~ 1 Mm, but AT is about 500 K at maximum, e.g., 
Fig.0. If one assumes that the energy related to the temperature 
fluctuations is transported vertically by (acou stic) w aves moving 
with a speed v ac of about 6kms _1 (cf. Sect. 13.5.51 later on), the 
generic chromospheri c energy requirement of about 4 kW m 4 of 
IVernazza et al.1 (Il976l) corresponds to about 0.67 J m~ 3 (horizon- 
tal dotted line in the top panel of FigO. The atmospheric en- 
ergy content from all three approaches (velocity, intensity, LTE 
inversion) intersects with this requirement at a height of about 
750 km, but drops rapidly below this value at higher atmospheric 
layers. The ratio of the generic energy flux of 4kWm~ 2 and 
the energy density multiplied with a constant sound speed of 
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6kms _1 is about ten at z ~ 1000 km and increases rapidly to a 
factor of several hundreds at layers above. 

The rapid drop of the energy density with height is caused 
by the exponentially decreasing gas density p gas that is over- 
plotted as well in the top panel of Fig. The increase in rms 
fluctuations with height is not able to compensate the decrease 
in p gas . We used two different gas densities for test purposes, on 
the one hand the tabulated values of the HSRA, Pg as , and on the 
other hand the gas density p in the modified HSRA model that 
uses a linear extrapolation above log r = -4. The differences 
between these two gas densities, however, do not significantly 
change the resulting energy density (compare the black solid and 
black dashed curve in the top panel of Fig. [6]). There is, however, 
a height range from about z = 100 km to 250 km, where the en- 
ergy of the temperature fluctuations from the LTE inversion stays 
roughly constant while the density reduces by a factor of two in 
the same height interval, similar to the height range of constant 
energy found for individual bright grains (see Fig. [TOllater on). 



3.4.2. (Radiative) energy losses 

The derivative of the energy content, again multiplied with v ac , 
gives the energy losses with height (middle panel of Fig. O. 
The region of roughly constant energy density can again be 
identified at about z ~ 200 km. Integrating the area of the en- 
ergy gradient for all heights above 750 km that can be consid- 
ered to contribute to the radiative losses in the Ca line core, 
i.e., the "chromosphere", yields a value of about 3.5 kWm -2 , 
which would match the generic chromospheric energy require- 
ment. This would, however, imply that no energy is left anymore 
for heating the layers above such as the corona since this amount 
of energy is assumed to be radiated away in all of the chromo- 
spheric spectral lines. 

For comparison with the gradient of energy losses 
with height, we used the value of chromospheric radia- 
tive losses derived f r om a set of numerical simulations by 
ICarlsson & Leenaartsl(l2012l CL12). Their Figs. 15 and 16 yield 
an average cooling rate of about 100 Wg" 1 in the chromosphere 
that can be converted to losses in energy density by a multipli- 
cation with the gas density as a function of height. The resulting 
energy gradient with height using a constant value of 100 W g" 1 
(red line in the middle panel of Fig. [6]) is larger than our esti- 
mate (black line in the middle panel of Fig. by a factor of 
about 2-3 between z ~ 200 - 300 km, and nearly identical above 
Z ~ 800 km. This difference might be caused to some extent by 
the fact that we only considered rms fluctuations of temperature 
in the energy estimate, whereas the value of CL12 refers to all 
radiative losses, i.e., losses related to the modulus of tempera- 
ture itself, not its variation. For a final comparison with CL12, 
we determined the energy loss rate per mass by dividing our 
energy gradient with height by the mass density (bottom panel 
in Fig. compare with Fig. 16 of CL12). The global trend of 
an increase in radiative losses from the photosphere towards the 
chromosphere is identical in both estimates with the prediction 
of a maximum at about 1 Mm. Less pronounced local maxima at 
z ~ 800 and 1300 km appear in both estimates, albeit in the case 
of our estimate their significance is doubtful. In total, a roughly 
constant value of radiative losses of about 100-200Wg _1 in 
the lower chromosphere fits to both our estimate and the one of 
CL12. 
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Fig. 7. Residual of an energy input of 100% at a height of 
240 km (solid black line). Top: linear scale. Bottom: logarithmic 
scale. The red pluses and blue diamonds give the corresponding 
values of FA 12 for a high (case a) and low (case b) magnetic fill 
factor. The green asterisks give the values of ST08 for acoustic 
waves (case a), the orange triangles thos e of gravity waves ( case 
b). The purple squares were taken from lFelipe et al.l (1201 ll) and 
pertain to magneto-acoustic waves in a sunspot umbra. 



3.4.3. Residual energy content 

Figure [7] shows the residual energy left at a given height when 
some amount of energy is dumped into the atmosphere at a 
height of z = 240 km. To clearly visualize the effect of the en- 
ergy losses with height, we provide both a linear and a loga- 
rithmic plot. The black curve corresponds to our result for the 
height variation of the internal energy density shown in the up- 
per left panel of Fig. but normalized to the energy density 
at z = 240 km. With the strong gradient of energy with height, 
only 5 % of the initial energy are left at a height of 500 km, sim- 
ilar to the reductio n of wave energy by 80 % at 400 km found in 
ICuntz et all (120071) . For comparison, we also overplotted the rel- 
ative energy fl uxes in the numerica l simulations of acoustic wave 
propagation of lFawzy et al.l (120121 FA12, first and last column of 
their Table 3). Their wave energy at 600 km was normalized to 
the corresponding fraction of our curve at that height. We also 
read off the value of the energy fluxe s for acoustic waves and 
gravity waves from lStraus et all (12008, ST08; their Fig. 3), both 
normalized to 100 % at z = 250 km. As last reference, the acous- 
tic energy flux in a wave simulation inside a sunspot umbra from 
iFelipe et all (120111 their Fig. 16) is also overplotted. 

A comparison of the various curves shows that our result 
on the variation of the energy density with height, which corre- 
sponds to the rms fluctuations of brightness temperature in the 
QS, matches best with all those derived for acoustic waves in an 
only weakly magnetized plasma (case a of ST08 and case b of 
FA 12). This implies that any gravity waves present in the atmo- 
sphere should have no signature in intensity spectra, and hence 
temperature, otherwise the inversion would have had to retrieve 



a significantly different slope of the energy density with height 
as given by case b of ST08. 

3.5. Individual bright grains in the quiet Sun 

Bright grains (BGs, cf. iRutten & Uitenbroeklll99ll) are one of 
the most prominent features in observations of the chromo- 
spheric Can lines, regardless whether two-dimensional (2D) 
imaging data or spectroscopic data are used. They are transient, 
strong intensity brightenings in chromospheric layers that ap- 
pear with a random spatial pattern, but a characteristic period 
of about three minutes. They are clearly related to some sort of 
energy transpor t towards and energy in crease in the solar chro- 
mosphere (e.g. JCarlsson & Steinlll997l CS97 in the following; 
BE08). CS97 studied the purely hydrodynamical case of BGs in 
numerical simulations and were able to show the photospheric 
origin of these features. 

Following the general approach used throughout our series 
of papers, we will treat all BG events in our data in the same 
way without trying to distinguish between a possible acoustic or 
magnetic origin, which would be possible using the spectropo- 
larimetric data in the photospheric Fe i lines at 630 nm. The si- 
multaneous polarisation measurements (cf. BE08, their Fig. 4) 
showed a stable magnetic environment wit hout clear indications 
of neither magnetic flux emerge nce (e.g.. iGomorv et aD 1201 Oh 
nor flux cancellation events (e.g jBellot Rubio & Beckll2005l) . 

To estimate the energy transport related to BGs, we stud- 
ied the occurrence rate, spatial extension, and energy content of 
individual BGs for an extrapolation of their total energy contri- 
bution to the chromospheric energy balance. We used the data 
of the time series taken at the disc centre (observation No. 2 in 
BE 12) for that purpose. We define BGs as any temporally iso- 
lated temperature enhancement in the lower solar atmosphere, 
i.e., also below chromospheric layers. 



3.5.1 . Temporal frequency 

To derive the characteristic occurrence frequency of BGs, we 
analyzed the temporal evolution of the temperature on the exis- 
tence of local extrema, i.e., local temperature maxima. The top 
panel of Fig. [8] shows the temporal evolution of the temperature 
excursions around the mean value at logr = -3.6 for one pixel 
along the slit as an example. The temporal evolution of the tem- 
perature is similar to the observed intensity at A = 396.83 nm, 
which roughly forms at that optical depth (red dash-dotted line). 
We then applied a method developed for finding polarisation 
lobes or local inten sity reversals in line profiles (e.g., Beck 2006; 
iRezaei et~ai]|2008l) to all pixels along the slit and all layers in 
log r to find the locations of local maxima in the temperature ex- 
cursions. To be counted as a local maxima, the value at a given 
instant of time has to be larger than the values on the two time 
steps before and after, i.e., it has to be the largest value in a range 
of ±40 sees. With this approach, the influence of noise on the de- 
termination is reduced. Given the typical duration and cadence 
of BGs (cf. Figs.[T0land[T3llater on), ±40 sees is a characteristic 
time span in which usually only a single BG happens. The tem- 
perature maps were filtered beforehand for all variations with pe- 
riods above 7.5 mins to suppress the long-term variation caused 
by the granular evolution. We implicitly identify in the following 
local temperature maxima with BGs. 

Along the slit, the number of maxima was comparable for all 
pixels, regardless if they were related to photospheric magnetic 
fields or not (middle panel of Fig. [8]). We find only a weak depen- 
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Fig. 9. Temporal evolution of the H2V intensity in a 2D FOV. The 
white vertical lines denote the area of one scan of about 1.5 Mm 
extent in x taken in about 21 sees each. The short white bar at 
(x, y) « (1 1, 6) Mm indicates a length of 1 Mm. 



sphere (cf. the top panel of Fig.[8]at t ~ 40 min). We repeated the 
determination of locations of local temperature maxima with a 
width of ±4 time steps (= ±80 sees) or after applying a running 
mean of five time steps width, which both yielded only about 
15 maxima on average. The dependence of the number of max- 
ima on optical depth did not change by the additional require- 
ment/smoothing. In case of analysing the intensity variation with 
time for local maxima in a ±40 sees range (top panel of Fig. [8]), 
three additional local maxima are found in comparison to the 
temperature curve. Thus, the most probable range for the number 
of temporally local intensity or temperature maxima, and hence 
BGs is from 15 to 22 in 3150 sees (0.005-0.007 sec" 1 ). 

The area fraction of BGs, i.e., the number of pixels along 
the slit that exhibit a temperature maximum at a given moment 
of time relative to the total number of pixels, is about 10 %. A 
similar estimate can be made from Fig. [9] that shows about one 
BG of 1 Mm diameter in each 2D FOV of about 1 Mm x 10 Mm. 



Fig. 8. Temporal occurrence frequency of BG events. Top: tem- 
poral evolution of temperature at logr = -3.6 for one pixel 
along the slit. The red dash-dotted line shows the intensity at 
X - 396.83 nm in arbitrary units for comparison. The vertical 
dashed lines indicate local maxima. Middle: temporal average 
of number of maxima along the slit. The dash-dotted red line 
shows the unsigned integrated Stokes V signal in arbitrary units. 
Bottom: average number of local maxima during 3150 sees as a 
function of log r. 



dence of the total number of maxima in the 3 150 sees duration of 
the time series on log r, with a slight increase in the occurrence 
rate at upper atmospheric layers (bottom panel of Fig. [8]). On 
average, we find about 20 maxima in temperature excursions at 
all slit positions and layers of log r , which yields a characteristi c 
cadence of about 160 seconds (cf. Ivon Uexkull & Kneedfl995l) . 
The lack of variation with optical depth may seem puzzling at 
a first glance, given that usually solar oscillations are separated 
into photospheric 5-min and chromospheric 3-min oscillations, 
but ascribing any period to not necessarily repeated events can 
also be misleading. For instance, the three BGs traced in Figs.fTOl 
and[T3llater on stay three events throughout all optical depth la y- 
ers (see also Kalkofen et al. 2010; Stangalini & Morettill20T2h . 

The number of 20 events in 3150 sees yields an occurrence 
rate of about 0.006 sec" 1 , wh ich is comparable to the values 
given in lTritschler et al.l (120071 their Fig. 7, 0.0075 sec -1 at max- 
imum). We note that our rate is an upper limit for the occur- 
rence of BGs because the simple identification method by local 
temperature maxima presumably has also counted some peaks 
that could result rather from observational effects such as resid- 
ual image jitter than from physical processes in the solar atmo- 



3.5.2. Horizontal spatial extent 

To estimate the spatial extent of individual BGs, we used the fact 
that the time series on disc centre (observation No. 2 in BE 12) 
actually consists of 150 repeated scans of 4 steps with 075 step 
width each, even if only the single step of each scan that was co- 
spatial in 396 nm and 630 nm was investigated in detail. Figure 
[9] shows the temporal evolution of the H2V intensity in a 2D 
FOV of 1.45 Mm x 14. 8 Mm over 651 seconds (=31 repetitions 
of scanning the same area). The three consecutive brightenings 
of, e.g., Fig. [K)|are located at (x,y) « (0 - 32, 6) Mm. Several 
BGs can be identified throughout the FOV, but their full extent 
in most cases is only covered along the slit in y whereas they 
are truncated in x when they were not almost centred inside the 
scanned area of 2" extent (e.g., the BG in the fourth repetition at 
y ~6Mm in Fig. [9]). 

From the section of the FOV shown in Fig. |3 a roundish 
shape with an average diameter of about 1 Mm can be attributed 
to BGs. To improve the statistics, we manually identified about 
100 BGs in the full FOV of all the time series. The BGs were 
required to be located outside of regions with strong polarisation 
signal and to be visible in at least three repetitions of scanning 
the same area. We then selected the BGs at the time of their 
largest brightness and took intensity cuts along the slit through 
their centres. We fit Gaussians to the intensity cuts to obtain the 
diameters of the BGs. The diameters ranged between 0.50 and 
1.83 Mm, with an average diameter of 1.17 Mm. A set of about 
30 BGs that happened on locations with significant polarisation 
signal showed a slightly larger average diameter of 1.36 Mm. In 
the following we will thus assume a typical radius of 0.5 Mm for 
a single BG, which agrees well with the characteri stic diameter 
of BGs of 1 -2" found in Can imaging d ata (e.g. JWoger et al.l 
l2006UWogerll2007UTritschler et al.ll2007l) . 
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Fig. 10. Estimate of the energy content of individual BG events. 
Top panel: temperature enhancements (left column) and internal 
energy enhancements (right column) for a pixel with BG activity 
(bottom row) and a quiet reference pixel (top). Bottom panel: av- 
erage internal energy density for the active (black) and the quiet 
pixel (red) vs. geometrical height. The blue line shows their dif- 
ference. 

3.5.3. Energy content 

To determine the energy content of individual BGs, one needs 
to derive the temperature enhancement related to them. One ad- 
ditional problem here is to separate the temperature contribu- 
tion of BGs from, e.g., the background (brightness) temperature 
that results from an atmosphere that is in radiative equilibrium 
with the incoming photospheric radiation. This background ra- 
diation also needs not be constant throughout the FOV because 
photospheric magnetic fields or variations in the granular evo- 
lution can modify it. We therefore first determined a minimum- 
temperature stratification T m { n (r,y) for each pixel along the slit 
defined by the lowest temperature at each layer of log r attained 
on the pixel y during the one-hour time series. A subtraction of 
this minimal temperature from the temperature stratification on 
the pixel y at a given moment of time then yields only positive 
temperature fluctuations. These fluctuations represent the mini- 
mum amount of additional energy above the lowest energy state 
that the atmosphere exhibited during nearly one hour, without 
excluding that already the lowermost energetic state was ener- 
getically above radiative equilibrium. 

To isolate the BG contribution to the energy and to remove 
any additional contribution from the background radiation, we 
first selected two different locations for a comparison, i.e., one 
pixel with a significant activity in form of a sequence of three 
prominent and well-defined BGs that appear consecutively, and 
one pixel that showed no clear BG or other activity (marked in 
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Fig. 11. Spatio-temporal occurrence of temperature maxima. 
Bottom row, left to right: intensity maps at A = 396.724, 396.822, 
and 396.851 nm. The horizontal white/black bars mark the loca- 
tions of the active/reference pixel used before. Middle row, left 
to right: energy density at heights of log r = -1.6, -3.3, and -4.6 
at the time of temperature maxima. Top: complementary map of 
energy density on all other locations. 



Fig- [IB- The top panel of Fig. [TCI shows the temperature fluctua- 
tions relative to r m i n (r, y) for the two pixels over a period of 65 1 
seconds, with the "active" pixel in the bottom row and the "quiet" 
reference pixel in the top row. The left column shows the tem- 
perature fluctuations AT and the right column the correspond- 
ing variations of the internal energy density AE that were cal- 
culated using Eq. (Q]). For the conversion to geometrical height, 
we assumed the tabulated relation between log r and geometri- 
cal height z of the HRSA to be valid as first order estimate also 
for individual spatially resolved temperature stratifications. 

To determine the characteristic energy content of a BG, we 
averaged the value of AE for the two pixels over 65 1 sees (bot- 
tom panel of Fig. [TO]). The difference of the two average values 
of AE reaches a maximum of 22.6 J m~ 3 at a height of about 
250 km. There are several arguments to use the difference at 
about this height as the actual energy input of one BG. Around 
this height, the average energy density is nearly constant for the 
active pixel, whereas it drops steeply at higher atmosphere lay- 
ers, indicating that the energy is inserted at low altitudes in the 
atmosphere. This is also evident from the temporal evolution of 
the temperature fluctuations, where the precursors of the BG ap- 
pear first low in the atmosphere in each case and then propa- 
gate upwards. The temporal evolution of the temperature (upper 
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Fig. 12. Temperature stratifications vs. geometrical height at the 
onset of two BGs corresponding to the vertical white dotted lines 
in the top panel of Fig.[l0j The horizontal dotted lines denote the 
FWHM of the BG events. The blue line shows a sinus curve with 
a wavelength of 400 km for comparison. 



panel of Fig. [TO]) also shows clearly that three BGs happen in- 
side the time used for the average, which are caused by three 
short-lived (< 100 sees) events of energy input into the low at- 
mosphere. The temporal fill factor - as estimated from the time 
intervals when the BGs can be clearly identified in the tempera- 
ture variation in Fig. [TO]- of all three BGs together in the time 
span considered is about 50%. The energy enhancement of a 
single BG thus should correspond to about twice the temporally 
average value, i.e., about 45 Jm~ 3 . 

To improve the statistical base of the result and to ensure that 
it is not (strongly) biased by exactly the train of BGs and the ref- 
erence pixel that we selected, we also calculated the statistics of 
BGs relative to non-BGs for the full FOV. To that extent, we used 
the mas k of lo cations with local temperature maxima determined 
in Sect. l3.5.ll The internal energy density of all locations, where 
local temperature maxima at three different optical depths oc- 
curred, is shown in the middle panel of Fig.[TT] whereas the top 
panel shows the complementary map of the energy density. The 
bottom panel shows the intensity at the wavelengths that were at- 
tributed to these optical depths. A comparison of the middle and 
upper panel reveals that the BGs with their apparent diameter of 
1 Mm are usually only a part of a wave front that extends over 
up to ten Mm (cf. lvon Uexkull & Kneel 19951) and induces tem- 
perature maxima along its full length. The temperature maxima 
mark the spines of intensity brightenings, with a similar rise and 
decay time of the BG events (top panel of Fig.Qj]). The temporal 
extent appears strongly compressed in the spatio-temporal map 
in comparison to the spatial extent. 

To determine the maximal energy enhancement caused by a 
BG, we calculated the average value of the energy density for 
all temperature maxima and for the complementary sample, i.e., 
all points in time and space that did not correspond to a local 
temperature maximum. For the optical depth of logr = -1.2 
(= 148 km), BGs have an average energy density of 135 J m~ 3 , 
whereas the corresponding value at logr = -2.2 (= 244km) is 
124 J m~ 3 . For the non-BG sample, the average energy density is 
89 J m" 3 and 80 J m~ 3 , respectively. We excluded not only the lo- 
cation of the temperature maximum, but also always the directly 
preceding and following time step in the calculation of the av- 
erage non-BG energy. The energy difference between BGs and 
non-BGs then is 46 J m" 3 and 44 J m" 3 at 150 km and 250 km 
height, respectively, matching the estimate above from individ- 
ual BGs. We will use a characteristic enhancement of the internal 
energy density of 45 J m~ 3 for one BG in the following. 
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Fig. 13. Vertical propagation speed of temperature perturbations. 
The background image shows the temperature variation of the 
active pixel of Fig. [TO] vs. geometrical height. The three dotted 
lines trace the propagation of the perturbations with time. 

3.5.4. Vertical extent 

To determine the vertical extent of BGs, in principle the curve 
of the energy difference between the quiet and the active pixel 
could be used (lower panel of Fig. [TO]), but we decided to use the 
temperature stratifications instead because of the larger ampli- 
tude of the temperature fluctuations. We converted the temper- 
ature stratifications at the onset of two BGs (indicated by the 
white vertical dotted lines in the lower left panel of Fig. [TO]) 
from optical depth to geometrical height (Fig. [12]). The BGs 
at this stage of their evolution show up as Gaussian-like tem- 
perature bumps centred at about z ~ 300 km. The full- width 
at half-maximum (FWHM) of the temperature enhancements is 
about 200 km in height, which would correspond to half a pe- 
riod of a sinusoidal wave with 400 km wavelength (blue line in 
Fig. [12]). Identifying BGs with acoustic waves with a period of 
3 or 5 min and assuming propagation with sound speed (about 
6 kms" 1 ) would predict significantly larger wavelengths of 1 and 
1.8 Mm instead. The propagation of BGs in the atmosphere is not 
purely vertical, but has also a horizontal component (cf. Figs. [9] 
and Fig. [TT] some of the BGs move along the slit with time), 
which modifies the apparent vertical wavelength. The angle be- 
tween local vertical and the wave propagation directions seems, 
however, to be small in most cases. 

With the information on the vertical extent and the effective 
enhancement of the internal energy density, one can derive an 
estimate of the total energy of a BG by E = V • AE. For a single 
pixel of 07 5 x07292 extent, assuming a vertical extent of 200 km, 
and using AE = 45 Jm~ 3 , one obtains a total energy of the BG 
of 6.91 x 10 17 J. Assuming a roundish structure with a radius of 
0.5 Mm, the total energy of a BG is about 7.1 x 10 18 J. 

3.5.5. Vertical propagation speed 

For estimating the vertical energy flux, the characteristic ver- 
tical propagation speed of BGs is needed. The LTE inversion 
provides the necessary information in form of temperature vari- 
ations in optical depth. We used again the tabulated relation be- 
tween r and z of the HSRA model to convert from optical depth 
to a geometrical height scale. Figure [13] shows the resulting tem- 
perature fluctuations at the active pixel vs. geometrical height. 
The precursors of the BGs leave inclined traces in the lower part 
of the atmosphere up to about 600 km height, whereas above the 
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temperature enhancements are nearly vertical, i.e., seemingly af- 
fecting all higher layers at the same time. The latter presumably 
does not indicate that the propagation speed increases strongly, 
but rather that the relation between r and z is not valid inside 
a shock front. CS97 concluded from their ID NLTE simula- 
tions that the emission in Can H at maximal brightness of a 
BG stems from a height range of about z ~ 1 to 1.4 Mm, where 
the strongest contributions to the emergent intensity originate in 
shocks. For the lower part and the onset of the BGs, the propaga- 
tion speed of the temperature perturbations can, however, well be 
measured (inclined dotted lines in Fig. [13]). The corresponding 
slopes denote vertical propagation velocities of about 5.0, 6.3, 
and 6.7 km s" 1 , respectively, c omparable to the sound speed (see 
also von Uexkull et al. I fl983h . We note that an estimate of the 
propagation speed and direction in all three spatial dimensions 
can be derived from an inversion of all steps of each repeated 
scan, dropping the restriction of using only the step that was co- 
spatial in the Ca n H and 630 nm spectra. 

The magnified view of Fig. [13] on the BG events allows 
one to estimate an additional characteristic, the typical duration. 
Above a height of 600 km, none of the BG events lasts longer 
than about 60 sees from its onset to its subsequent complete dis- 
appearance. 

3.5.6. Energy flux related to BGs 

The temporal evolution of the intensity or temperature shows 
that the energy of BGs is removed by both propagation and ra- 
diation in a short time. There are, however, different ways to 
estimate the vertical energy flux per unit area. Using the energy 
content of a single BG of 6.91 x 10 17 J per pixel of 075 x07292 
as estimated above and assuming that all this energy disappears 
again before the onset of the next BG 160 sees later yields an 
energy flux of 56.3 kWm -2 , but without specifying into which 
spatial direction (or into which degree of freedom) the energy 
is dissipated or disappears. The average energy enhancement of 
22.6 J m" 3 of the active relative to the quiet pixel (Fig.[l0]) yields 
an energy flux of 41.4kWm~ 2 from both the influx and dissi- 
pation of the energy of three BGs in a time span of 651 sees. 
Excluding a downward energy transport and otherwise equal 
probability for an energy transport to the remaining five possi- 
ble spatial directions, the upwards vertical energy flux would be 
about (41.4-56.3)/5 ~ (8.3- 11.3)kWnr 2 . This estimate uses 
the characteristic temporal evolution in the derivation, i.e., about 
160 sees cadence between BGs or three events in 651 seconds, 
respectively, and refers to the energy flux between about 150 to 
350 km height. 

However, even if individual BGs induce a vertical energy 
flux of up to a few tens of kW m~ 2 in the low solar atmosphere 
(z < 400 km), it seems that at z ~ 1 Mm far less than 5 % of 
the energy flux are left over (Fig. [7]>. We remark that in our case 
the underlying curve that predicts this behaviour is directly pro- 
portional to the temperature rms with height and the gas density 
(Eq. []]). Modifications in the stratification of either of the two 
quantities therefore will change the residual energy as well. 



4. Discussion 

The reliability of our analysis results depends to some extent on 
the LTE assumption taken. We therefore first return to a discus- 
sion of observed and theoretical profiles that is partly indepen- 
dent of the (N)LTE question. 
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Fig. 14. Comparison of average observed and best-fit profiles 
with theoretical NLTE profiles. Bottom panel: spectra from a 
time series in the QS on disc centre. Top panel: spectra from 
an AR map off centre. Thick black: average observed QS pro- 
file. Thin purple: average best-fit QS profile. Blue: LTE profile 
of modified HSRA. Red/orange: NLTE FAL-C/A profile. For the 
AR, average profiles of a plage area (blue dotted) and an umbra 
are over-plotted (orange dash-dotted). The vertical dotted line 
denotes the wavelength of intensity normalization. 

4.1. Comparison of average profiles and corresponding 
temperature stratifications 

Figure [14] shows the average observed and best-fit profile in QS 
for our time series and for one of the AR maps, an average um- 
bral and plage profile from the same AR map, the profile synthe- 
sized in LTE from the modified HSRA model without a tempera- 
ture rise, and two profiles synthesized i n NLTE from the FA L-A 
and FAL-C te mperature stratificat ions (iFontenla et al.ll2006l) us- 
ing the code of lUitenbroekl (1200 ll) . The spectral sampling of the 
FAL NLTE profiles in the synthesis was 0.75 pm without any ad- 
ditional spectral degradation (cf. the next section). To facilitate 
an easier comparison of the variations in the very line core in 
this figure, all profiles from the AR map were normalized to the 
same intensity at 396.9 nm. 

Both NLTE FAL profiles exceed the observed emission 
peaks in the QS on disc centre (lower panel of Fig. [14]). The 
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Fig. 15. Comparison of temperature stratifications. Bottom'. QS 
maps. Top: one of the AR maps. Thick black: average QS tem- 
perature in the LTE fit. Blue: modified HSRA model. Dashed: 
original HSRA model. Red: FAL-C. Orange: FAL-A. For the 
AR, the average temperature stratification of a plage area (blue 
dotted) and of the umbra are over-plotted (orange) instead. 



FAL-C profile exceeds even the emission in the QS region of the 
AR map (top panel). To exclude a possible influence of the large- 
scale averaging that could reduce the amplitude of the emis- 
sion peaks through the presence of spatially/temporally varying 
Doppler shifts, we calculated a QS profile that was averaged 
only over 112 spectra, i.e., one slit position in the time series. 
This small-scale average QS profile showed a slightly increased 
amplitude of the H2V emission peak, but the change stayed mi- 
nor and the intensity still fell far short of that in the NLTE FAL 
profiles. Both the observed and average best-fit profile, however, 
clearly exceed the profile corresponding to the modified HSRA 
atmosphere. Average profiles in the QS therefore fall between an 
profile from an atmosphere in radiative equilibrium and an atmo- 
sphere with a chromospheric temperature rise such as in FAL-C 
or FAL-A. For the case of the AR map, the profiles change sig- 
nificantly. The observed average QS profile has two clear emis- 
sion peaks, and the average best-fit profile in the QS starts to 
resemble the NLTE FAL profiles. We remark that inside the AR 
map, the "QS" regions contain far more magnetic flux than found 
in the QS at disc centre, and furthermore that the AR observa- 
tions were taken at a heliocentric angle of about 50°. The emis- 
sion in average plage and umbra exceeds the FAL NLTE profiles 
by a factor of about two. 

The corresponding temperature stratifications in Fig. [15] re- 
flect closely the differences in the profiles themselves. The mod- 
ified HSRA lies below the average temperature of the QS on disc 
centre retrieved by the inversion, and the two FAL models above 
it, whereas the stratifications obtained by the inversion for the 
AR region plage and umbra follow to some extent the chromo- 
spheric temperature rise in the FAL models or the original HSRA 
model (top panel in Fig.[T5l). 

The LTE assumption for sure has a tremendous effect on 
the retrieved temperature stratifications, but the main prob- 
lem consists already in reconciling the observed average spec- 
tra in QS with those predicted by the NLTE FAL profiles 



(lower panel of Fig. [14]). Appendix lAl discusses possible spec- 
tral broadening mechanisms that could reduce the amplitude of 
the emission peaks in observed spectra. We find that any rea- 
sonable instrumental broadening, or broadening by unresolved 
micro- and macroturbulent velocities with reasonable velocities 
(< lOkms -1 ) are insufficient to reduce the emission peaks from 
the NLTE FAL models enough to fully match the observed spec- 
tra, leaving as o nly option a clear re duction of the chromospheric 
temperature (cf. lRezaei et al.ll2008l) . 

4.2. Difference between temperature stratifications in quiet 
Sun and active regions 

Even if it is well known that the emission in the Can H 
line core inc reases in the presence of ( photospheric) magnetic 
fields (e.g., [Simon & Leightonl 1 1964 



Linskv & Avrett 1970 



| Liu & Smithlll972l:lDunn & Zirkedll973[lMattig & Kneed 1978. 
ISchrijver et al. Rezaei et al H2007IK the difference between 
temperatures in QS and AR retrieved by the inversion is surpris- 
ingly large. The LTE approach used in the inversion retrieves 
temperature stratifications with a clear chromospheric rise in the 
AR, even though it supposedly underestimates the temperature. 
That the underestimation comes with a caveat is caused by a di- 
chotomy in the (N)LTE problem. In the LTE assumption, (gas) 
temperature and intensity are directly coupled. In case that an 
increase in emission is not caused by an increase in the gas tem- 
perature but by NLTE effects on the emitted radiation field, the 
LTE inversion will attribute it still to a temperature change and 
will yield a too high gas temperature. In case that the emission is 
caused by an increase in the gas temperature, the LTE approach 
will yield a too low temperature, because it assumes that the cor- 
responding energy will be evenly distributed over all degrees of 
freedom, whereas in NLTE a significant fraction can go into ad- 
ditional ionization or other energy sinks without increasing the 
emitted radiation. 

For active regions, the temperature structure in the solar 
atmosphere seems to be more stable than in the QS. While, 
e.g., umbral profiles show clear indications for upwards prop- 
agating wave patterns produ cing shock waves in the chromo- 
sphere like in the QS (e.g., iLitesI 119861: ICenteno et all 120061: 
iFelipe et al J 12010k Can H spectra in umbrae never loose the ir 
emission peaks (Mattig & Kneed 1 19781: Kneer et all U981blS 
iRouppe van der Voort et al. Felipe et al. r i2oToi) . m con- 
trast to QS profiles that oscil late betwe e n low or no emission 
and intensity reve rsals (e.g, iLiul 11974 ICram & Darnel fl983l 
iRezaei et al J20081 BE 12). The strong v ariability in QS was also 
found in numerical si mulations (CS97, We demeyer et al.l l2004: 
iLeenaarts et al.ll20TIk but the corresponding prediction for ac- 
tive regions or the interior of sunspots is still pending. Even if 
we cannot provide any conclusive proofs, the results of our LTE 
inversion support a more stable gas temperature structure in ac- 
tive regions with a clear chromospheric temperature rise. This 
does not exclude that the line shape is caused by NLTE effects, 
but suggests that the underlying temperature stratification seems 
to be less influenced by temporal variations than in the QS. 

The modulus of the difference between QS and active re- 
gions even in our LTE approach indicates that heating processes 
related to strong magnetic fields are up to five times stronger 
than those related to acoustic waves in the QS: the (spatial) tem- 
perature rms in the AR is about five times larger than the spa- 
tial or temporal rms in QS (Fig. [2j and the enhancement of the 
QS temperature stratification above the modified HSRA model 
(AT < 1000 K, Fig. [15]) is roughly five times smaller than for 
the average plage (AT ~ 4000 K). Because in most of the stud- 
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ies of low chromo spheric lines such as Ca n H the presence of 
magnetic fields was neglected, their inclusion could be the turn- 
ing point to reconcile the high variability of, e.g., Can H with 
the low variability of chromospheric lines observed by SUMER 
(iCarlsson et al.lll997l) . We point out that the latter authors also 
found a permanent presence of emission in the SUMER lines in 
QS, in conflict to Can H spectra in QS, but in full agreement 
with the behaviour of Ca n H spectra in the umbra. This would 
fit all together with the assumption that the SUMER lines form 
in a magnetized atmosphere (like Can H in the umbra) in QS, 
e.g., above a magnetic canopy with a strong damping of imping- 
ing acoust ic oscillation, wh ereas Ca n H forms below the canopy 
in QS (cf. lBeck et aDl2Q08L their Fig. 22). We note also that the 
same effect is suggested by the numerical simulations of FA 12. 
The average temperature stratifications shown in their Fig. 6 do 
not exhibit a clear temperature reversal below a height of 1 Mm, 
and additionally show a strong increase of the temperature when 
the magnetic fill factor is increased by an order of magnitude, 
caused by the related reduction of the gas density in a partially 
evacuated magnetized plasma. 

4.3. Energy content of bright grains and related energy flux 

We estimated the energy content of the precursors of chromo- 
spheric bright grains. We find a spatial extent of about 1 Mm 
by 0.2 Mm in height, with an increase of the internal energy by 
about 45 J m" 3 and a total energy of about 7 x 10 18 J per event. 
The temperature enhancements propagate upwards wi th about 
6km s" 1 , comparable to the local sound speed (see also IWorralll 
120121) . With the typical temporal cadence of 160 s for individ- 
ual events this yields a vertical energy flux of about lOkWm" 2 
at a height of about 250 km. Some insecurity in this calcula- 
tion is related to the characteristic cadence because it is sensi- 
tive to how a local maxima is identified and counted. Our de- 
termination yields from 15 to 22 maxima in 3150 sees, which 
gives roughly a 20 % error on the energy flux. Our specific data 
is less suited for the determination of the occurrence rate of BGs 
than time series of intrinsically two-dimensional data. The num- 
ber of about lOkWm" 2 falls in line with the energy fluxes be- 
tween 3 to 12kWm" 2 at ab out 250 km height determined from 
ve locity oscillations in, e .g., iBello Gonzalez et a D J20Q9L [2oTob 
or lMalherbe et al.l (1201 2l MA 12). The latter authors also deter- 
mined the total energy contained in what they call "acoustic 
events", i.e., a phase of upwards directed energy transport that 
can cover a few wave trains and can last up to lOmin (cf. Fig. 9 
of MA12). They found total energies of about 1.9 x 10 19 J per 
acoustic event, which fits to the example used in our Sect. 13.51 
with three BGs of in total about 2.1 x 10 19 J in 651 seconds. 

MA 12 found an occurence rate of acoustic events of about 
5.7 x 10" 16 s" 1 m" 2 . With 20 temperature extrema in 3150 sees 
on each pixel of 073 x 075, and a spatial fill fraction of BGs of 
about 10 % along the slit, we obtain a rate of 8.3 x 10" 15 s" 1 m" 2 , 
or about 0.5 mm" 1 Mm" 2 , respectively. The difference between 
the two rates is about a factor of 15, where our rate should 
be larger by a factor of 3-5 because of the existence of 
multiple temperature maxima during an acoustic event as de- 
fined by MA 12. The remaining difference could be caused 
by the lower spatial resolution of our data, where the area 
of BGs will be slightly overestimated because of the spatial 
smearing. The recent ly inaugurated 1.5-m GREGOR telescope 
(ISchmidt et al.l l2012h will offer the possibility to obtain high- 
resolution, high-cadence spectra of the Can IR lines around 
850 nm with the Gregor F abry-Perot Interferometer (GFPI, e.g., 
iPuschmann et alJ l2012allbllch and of Can H with the BLue 



Imaging Solar Spectrometer (BLISS, IPuschmann et al.l l2012cl) 
in the near future. 

Photospheric energy fluxes in the range of a few kWm" 2 
would suffice to replenish the generic chromospheric energy re- 
quirement of 4 kW m" 2 , if - and only if - this amount of energy 
would finally reach the chromosphere. The estimates of the en- 
ergy losses with height, however, predict that less than 5 % of 
an energy input at 250 km height reaches layers above 400 km. 
This suggests that the chromospheric radiative losses cannot be 
balanced by the energy of propagating waves or temperature per- 
turbations despite a sufficiently large reservoir of mechanica l en- 
ergy in the photosphere (but see also iPetrukhin et al]|2012l) . We 
note that the temperature stratifications and the rms temperature 
fluctuations retrieved by our inversion only predict on average a 
drastic loss of energy with height, but do not directly provide in- 
formation about where the energy is going to. To determine the 
exact fraction of the energy dissipated by radiation, thermal con- 
duction, or changes of the ionization fraction requires a detailed 
study of individual events. 

4.4. Visibility and extent of magnetic canopies in temperature 
maps 

Up to here, we discussed the properties of chromospheric bright- 
enings and the chromospheric thermal structure in the quiet Sun 
as a purely hydrodynamical phenomenon without considering 
the presence of magnetic fields. The simultaneous polarimet- 
ric observations of the 630 nm channel of POLIS, however, also 
provide information on polarisation signals, hence photospheric 
magnetic fields. We can thus investigate the influence of the 
magnetic topology on the thermal structure in the quiet Sun as 
well. 

The typical behaviour of the Can H line with its repeated 
change from reversal-free absorption to emission and vice versa 
in "field-free" regions, i.e., on locations that are polarisation- 
signal-free down to amplitudes of about 10" 3 of I c , and with 
a persistent stationary emission on locations with photospheric 
polarisation signal fits to the idea of magnetic field lines as a 
mediator between a stationary and transient emission pattern in 
Can H. For photospheric magnetic flux concentrations, geomet- 
rical models can be constructed by setting field strength and to- 
tal magnetic flux in the photosphere, and requiring (magneto- 
hydrostatic pressure balance with their surroundings. This 
yields the classical "expanding flux tubes" whose diameter in- 
creases with height in the atmosphere (e.g. JSolanki et al.ll 19911 
S091). Now, if the magnetic fields have a profound influence 
on the chromospheric temperature structure, they have to show 
up in the temperature stratifications retrieved by the inversion, 
regardless of the LTE - NLTE problem. 

Figure [16] shows that this is (clearly) the case. Panel a shows 
the temporally averaged temperature stratifications along the 
slit vs. optical depth. The also temporally averaged unsigned 
Stokes V signal is overplotted as a black line as a reference of 
photospheric magnetic fields that were roughly persistent over 
the duration of the time series of about one hour. Already in 
the temporally averaged temperature, flux-tube-like features are 
found co-spatial to the largest values of |V|. To increase the 
visibility, and because the main effect of the magnetic fields 
is an increase of temperature relative to field-free surround- 
ings, we subtracted the temporally and spatially averaged tem- 
perature stratification from the temporally averaged tempera- 
tures (panel b). In that case, temperature enhancements of +50, 
+ 100, and +150K over the average temperature almost exclu- 
sively occur on locations of photospheric magnetic fields. For 
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Fig. 16. Spatial temperature and velocity maps of the time series. Panel a: temporally averaged temperature <T> t vs. optical depth. 
Panel b: <T> t minus temporally and spatially averaged temperature stratification <T> vs. geometrical height. Panel c\ temperature 
at one scan step minus <T> vs. geometrical height. The z-axis in the upper two panels is stretched by a factor of five relative to the 
x-axis. The black lines at the lower border of each panel denote the total circular polarisation along the slit in arbitrary units. The 
blue lines show one flux tube model of S091. The orange/red/purple contours in panel b (c) mark temperature excesses of 50, 100, 
150 K (100, 200, 300 K). 



the two strong magnetic flux concentrations at x ~ 14.5 Mm 
and 35 Mm, the temperature enhancements at chromospheric 
layers above 500 km given by the +50 K contour form exactly 
the "wine-glass" shape predicted and used in many photospheric 
flux tube models (SQ9L iPizzo et al.lll993l: iGadun et al.l l200ll: 
IUitenbroekll20Tl iKhomenko & Colladosll2012l their Fig. 10). 
In both cases, the temperature enhancements yield a nearly hor- 
izontal (temperature) canopy at a height of about 620 km that 
extends betwe en 1 and 5 Mm horizontally away from the central 
axis (see also Ide la Cruz Rodriguez et aLll2013l) . The tempera- 
ture enhancements delineate vertical structures at the end of the 
canopies. 

The restricted horizontal extent of the temperature enhance- 
ments is better seen in panel c that shows the temperature strat- 
ifications for one step of the time series. Close to photospheric 
polarisation signals, the nearly horizontal canopies in tempera- 
ture appear again prominently, especially from x ~ 6 to 16 Mm 
where three neighboring flux tubes create a continuous temper- 
ature canopy, but in-between separated flux tubes areas of sev- 
eral Mm extent with a strongly reduced temperature show up 
(x ~ -7 Mm, x ~ 23 -33 Mm). In simil ar- sized areas, rever sal- 
free Can H spectra were found before (iRezaei et al.ll2008l) . In 
total, this implies that on the one hand magnetic canopies have 
a direct effect on the temperature stratification even in a single 
snapshot of the solar atmosphere, and on the other hand there are 
no clear indications for volume-filling magnetic fields in the for- 
mation height of Can H in the temperature map of Fig. [T6] The 
appearance changes a bit in the temperature map for the large - 
area scan at high layers (observation no. 1 of BE09, cf. Fig. 0). 



Figure [16] clearly demonstrates that photospheric magnetic 
fields are a second acting agent for the chromospheric tempera- 
ture structure in addition to the transient BGs. A more detailed 
study of the modification of the thermal structure at or close to 
the locations of photospheric magnetic fields is, however, be- 
yond the scope of the present paper. 

4.5. Applicability of LTE Inversion 

Because of the known problems of assuming LTE in an inver- 
sion of any chromospheric line, it might be worth to note a few 
examples for which purposes the LTE inversion still serves with- 
out any, or with only weak restrictions. We remark that mainly 
the modulus of the temperature derived from the LTE inversion 
has a high inaccuracy because of the missing instantaneous cou- 
pling between gas temperature and radiation field in NLTE con- 
ditions, and that the conversion from optical depth to geometri- 
cal height suffers from the lack of hydrostatical equilibrium. The 
latter could, however, be obtained post-facto with a presumably 
minor impact on the resulting synthetic spectra. 

The detectability of canopies in temperature maps - that are 
implicitly assumed to correspond to the theoretically expected 
magnetic canopies - does depend only slightly on the modulus 
of temperature. Our results clearly show that the signature of the 
canopies can be detected in the surroundings of photospheric 
magnetic flux concentrations up to a few Mm (cf. F ig.[T6l). This 
is in contrast to the results of lJudge et all (1201 Oh who found 
a strong "confinement" of intensity enhancements, and hence 
magnetic flux in the chromosphere. The mismatch of these two 
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results shows that a quantitative in-depth analysis of spectra, 
even with the limiting LTE assumptions, is crucial for tracing 
the topology at chromospheric layers. 

For estimating the energy of individual BGs, we used a 
height range, or equivalently, a wavelength region in the line 
wing, where NLTE effects can be assumed to play a minor role 
(z < 500 km). These results are thus rather robust and would 
presumably not change significantly when a different inversion 
method would be used. The advantage of the inversion of the 
Ca ii H spectra now is to provide on the one hand also a quanti- 
tative guess for the corresponding properties at the atmosphere 
layers above, which cannot be accessed by any photospheric 
line, and on the other hand to provide clear information on 
the temporal evolution, such as the vertical propagation speed 
(Fig. [13]), or the exact propagation direction from combining hor- 
izontal motions in the 2D sections of the FOV (Fig. [9]) with the 
vertical propagation speed. 



5. Summary and conclusions 

In a LTE inversion of Ca n H spectra, we cannot retrieve on aver- 
age a chromospheric temperature rise in quiet Sun (QS) regions, 
whereas in active regions (ARs) this is the case. The difference 
between QS and ARs can be directly traced back to the observed 
input profiles themselves, independent of the LTE assumption 
in the inversion. We confirm the findings of BE 12 of a layer 
at about 150 km height with an increased occurence of high- 
intensity, and hence high-temperature events near and in mag- 
netic features. A comparison of line-core and line- wing images, 
and the correlation coefficient between intensities and temper- 
ature vs. optical depth suggests that the increased skewness of 
the intensity histograms on magnetic locations is caused by a di- 
rect signature of magnetic fields in temperature at two different 
atmosphere layers, i.e., both in the chromosphere and the mid 
photosphere. The magnetic network is inconspicuous at the lay- 
ers in-between. The appearance of the network in the mid pho- 
tosphere could be caused by the "hot wall" effect in combination 
with the shif t of the optical depth scale inside magnetic fields 
(ISpruitll 19761) , whereas the appearance in the chromosphere in- 
dicates an energy deposit instead. 

The energy contained in the rms fluctuations of temperature 
and the large difference between QS and AR profiles and tem- 
perature stratifications indicate the importance of chromospheric 
heating processes related to the presence of magnetic fields. The 
physical processes that could yield such a heating should be in- 
vestigated in more detail in the future and included in the anal- 
ysis of chromospheric lines. The temperature fluctuations in the 
QS fall short of the requirements needed to maintain a stationary 
chromospheric temperature rise. 

The properties of individual (heating) events related to chro- 
mospheric bright grains (BGs) can be determined during all of 
their existence in the LTE inversion results, from the BG precur- 
sors in the photosphere to their culmination in the lower chromo- 
sphere. The mechanical energy content and the vertical energy 
flux in the photosphere would be sufficient to fulfill the generic 
chromospheric energy requirement of about 4kWm~ 2 , but only 
a fraction well below 5 % of the mechanical energy reaches the 
chromosphere, falling short by far of the generic requirement. 
The decrease of the internal energy density with height agrees 
with that expected for acoustic waves. 

The topology of magnetic fields in the chromosphere is 
seemingly directly reflected in the temperature structure at layers 
above a height of 500 km also in an LTE analysis, as indicated by 



the appearance of canopy structures around magnetic field con- 
centrations both in average and in instantaneous x-z temperature 
maps. With all of the shortcomings of the LTE approach, this 
quantitative analysis of large-scale data sets in a chromospheric 
spectral line provides a wealth of interesting features that can be 
used to address open questions in solar chromospheric physics 
of which this paper still has not covered all despite its extent. 
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Fig. A.l. Convolution of the FAL NLTE profiles with different 
spectral PSFs. Bottom: convolution of FAL-C with Gaussians of 
o~ = 3.4 pm (blue) and 9.4 pm (purple). The black line shows the 
observed average QS profile. Left column: line wing. Right col- 
umn: line core. Top: convolution of FAL- A with a wavelength- 
dependent Gaussian with cr between 3.4 and 9.4pm. 



Appendix A: Effects of spectral broadening 

The effective spectral resolution of the POLIS spectrograph can 
lead to a loss in amplitude of the emission peaks by spectral 
smearing because of their small spectral extent. The default 
spectral resolution of POLIS is 220.000@400nm (= 1.8pm), 
but because of the on-chip binning by a factor of two, the 
resol ution is sampling-lim ited to about 3.8pm. Convolving an 
F TS (iKurucz et al.lll984l) refe renc e spectrum with the method 
of lAUende Prieto et al.l (120041) and lCabrera Solana et al.l (120071) 
to match the observed average Ca profile yielded the need for 
a Gaussian with cr ~ 1.7 pm. We note that the FTS spectrum 
contains line broadening by macroturbulent velocities because 
of being a spatially averaged QS spectrum. 

For an estimate of the maximal spectral smearing between 
the observed average QS profile and the FAL-C NLTE profile - 
synthesized without a macroturbulent velocity - by instrumen- 
tal effects, we first convolved the FAL-C NLTE profile as ref- 
erence by a Gaussian to match the shape and width of the line 
blends in the Ca wing in the average observed spectrum. The 
lower left panel of Fig. IA. II shows that this requires a value of 
cr of about 3.4 pm, which is in rough agreement with the predic- 
tion from the sampling limit. A convolution of the full spectrum 
with this value, as an estimate of the instrumental and velocity 
broadening in the atmosphere, reduces, however, the amplitude 
of the emission peaks of the FAL-C NLTE profile only slightly 
(lower right panel of Fig. IA. lb . To force the emission peaks in 
the synthetic profile to a shape roughly resembling the observed 
profile requires a convolution of FAL-C with a Gaussian with 
cr = 9.4 pm. The line blends in the wing, however, clearly ex- 
clude such a large instrumental broadening because in the FAL- 
C NLTE profile convolved with cr = 9.4 pm (purple line) the 
blends are twice as broad and half as deep as in the observed 
spectrum (lower left panel of Fig. IA.ll) . 

We then decided to try a wavelength-dependent width of the 
Gaussian used as spectral point spread function (PSF). The width 
was set to 3.4pm in the line wing up to A = 396.75 nm, gradu- 
ally increasing to 9.4 pm at A = 396.85 nm, and reducing down 
to 3.4pm again at A = 396.94 nm. As it seemed impossible to 
match the FAL-C NLTE and the observed profile even with a 



large broadening in the very core, we applied this variable con- 
volution to the FAL- A NLTE profile instead. The temperature 
stratification of FAL- A has a weaker chromospheric temperature 
rise than FAL-C (Fig. [15]) and the corresponding NLTE profile 
differs significantly from FAL-C, especially in the residual line- 
core intensity (Fig. [14]). The approach with the variable spec- 
tral PSF maintains the match in the line wing (upper left panel 
of Fig. IA. lb and leads to an acceptable agreement between the 
emission peaks of the convolved FAL-A NLTE profile and the 
observed spectrum (upper right panel). 

The wavelength-dependent Gaussian spectral PSF would 
correspond to a microturbulent velocity v miC that changes from 
2.6 kms -1 in the photosphere to 7.1 kms -1 in the chromosphere. 
The NLTE synthesis already included a height-dependent micro- 
turbulent velocity of more than 5 kms" 1 for logr < -5.5. Thus, 
the effective microturbulent velocity necessary to match the av- 
eraged observed and the FAL-A NLTE profile should be nearly 
twice as large as the initial value of v miC , and has to be present 

already at lower layers of log r 3 to yield the correct shape 

of the Hiv and Hi R minima next to the emission peaks. We note 
that similar values of v miC would be also required to match the 
FTS atlas and the FAL-A NLTE profile, because the line shape 
in the FTS is similar to the average QS spectra shown here. If 
such a high value of v miC is reasonable for spectra at about 1" 
spatial resolution needs a further investigation. Additionally, the 
increase of v m ; c at comparably low layers of log r also already 
affects the Fei line blend at 396.93 nm strongly, broadening it 
significantly beyond its width in observed spectra. Doubling v m / c 
therefore seems not to be a valid option to match observed and 
theoretical NLTE p rofiles, leaving a re duction of temperature as 
the only option (cf. lRezaei et al.ll2008l) . 



